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We study charged boson stars in scalar-tensor (ST) gravitational theories. We analyse the weak 
field limit of the solutions and analytically show that there is a maximum charge to mass ratio for 
the bosons above which the weak field solutions are not stable. This charge limit can be greater 
than the GR limit for a wide class of ST theories. We numerically investigate strong field solutions 
in both the Brans Dicke and power law ST theories. We find that the charge limit decreases with 
increasing central boson density. We discuss the gravitational evolution of charged and uncharged 
boson stars in a cosmological setting and show how, at any point in its evolution, the physical 
properties of the star may be calculated by a rescaling of a solution whose asymptotic value of the 
scalar field is equal to its initial asymptotic value. We focus on evolution in which the particle 
number of the star is conserved and we find that the energy and central density of the star decreases 
as the cosmological time increases. We also analyse the appearance of the scalarization phenomenon 
recently discovered for neutron stars configurations and, finally, we give a short discussion on how 
making the correct choice of mass influences the argument over which conformal frame, the Einstein 
' frame or the Jordan frame, is physical. 
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I. INTRODUCTION 



^ ■ Boson stars are localised, asymptotically flat configurations of gravitationally bound zero temperature bosons. 
— \ Mathematically, the boson field is described by a complex wave function ijj whose Lagrangian possesses an internal C/(l) 
symmetry that gives rise to a conserved charge N , interpreted as the total number of bosons. The time dependence 
of Tp does not appear in the field equations, and the solutions describe bound cigenstates of ^ with a number of nodes 
that increases with the energy. The first boson star solutions were found by Kaup |Q| and independently by Ruffini 
I & Bonazzola Q|, who studied spherically symmetric stars in General Relativity (GR). They found that the solutions 
. were qualitatively similar to those describing neutron stars and white dwarfs, although they were of much smaller 

■ mass. Colpi, Shapiro & Wasserman extended their work by examining GR boson stars whose matter Lagrangian 
includes a quartic self-interaction term. This additional term contributes to the pressure of the star, increasing both 

^ its mass and particle number. More general self-interaction terms can be considered. The first research on this subject 
bJ[)' was undertaken by Lee and co-workers who, in a series of papers, studied the properties of non-topological solitons. 
' These are boson stars whose self-interaction term allows localised, non-singular solutions to exist even in the absence 
• ^ : of gravity. To accomplish this, the self-interaction must have attractive terms and must be at least sixth order in the 

■ boson field amplitude. The work of Lee et al. was performed independently of the main body of boson star research. 
^ 1 For the case of non-soliton stars, it is not immediately obvious how the inclusion of a non-quartic self-interaction 

■ " " ' affects the properties of the star. This problem will be addressed elsewhere |5[|; for our purposes, it will be sufficient 
to only consider the usual quartic self- interaction term. For reviews see Ref. 

Charged boson stars in GR were first studied by Jetzer & van der Bij 0] who found that the inclusion of the charge 
increases the star mass and particle number, in much the same way as the inclusion of the quartic self-interaction. 
In addition, they found that there was an upper limit on the charge to mass ratio of the bosons above which no 
non-singular solutions could be found: stars made up of bosons whose charge is beyond this limit generate a repulsive 
Coulomb field that overcomes the gravitational attraction, regardless of the number of bosons that make up the star. 
The stability properties of charged boson stars in GR were studied by Jetzer |H before the advent of the use of 
catastrophe theory in stellar equilibrium. Also, Jetzer and coworkers |^ have discussed the influence of charged boson 
stars on the stability of the vacuum and have considered their contribution to cosmological dark matter. 

Recently, there has been renewed interest in alternative theories of gravity, in particular scalar-tensor (ST) theories, 
in which the usual metric gravitational field is augmented by a scalar field (p which couples to the curvature via a 



*Email: A.W.Whinnett@qmw.ac.uk 

t Email: dtorres@venus .fisica. unlp .edu. ar 



1 



parameter uj{4>). The strength of this couphng increases with l/u. The simplest of these is Brans Dicke (BD) theory 
pof in which u is constant. These present studies of ST theories are motivated by the fact that they appear as the 
low energy limit of string theory (the simplest string effective action being the w = — 1 BD action) and that 
scalar gravitational fields arise from dimensional reduction of higher dimensional theories (see Ref. and citations 
therein). These theories may better describe gravity in the early Universe and, in addition, several models of inflation 
are driven by the same scalar field of ST gravity jlSj. Finally, ST theories provide a self-consistent framework in 
which one can study the possible variation of the gravitational coupling G. 

The strongest constraints on ST theories are usually assumed to come from the solar system weak field tests [ p4[ . 
As is well known, observations constrain the BD parameter to have a value of > 500 to within la. However, more 
general ST theories can satisfy the rather stringent constraints determined at the current epoch, while still differing 
considerably from GR in the past . In addition, it has been shown that strong non-perturbative ST gravitational 
effects may arise in ST neutron star configurations at the present epoch, even when the parameters of the theory 
accord well with observational tests . 

The current popularity of alternative models of gravity has prompted several authors to study boson stars in the 
framework of ST theory. The first of these were Gunderson & Jensen jl^ , who considered BD boson stars with cu = 6. 
More general ST boson stars were considered by Torres and later by Comer & Shinkai [|T9| . These studies showed 
that the inclusion of the field introduces no qualitative change in the solutions, although in all cases, boson stars 
have a slightly smaller mass than their GR counterparts. 

Spherically symmetric charge distributions in ST theory have also been investigated before. Singh & Usham |2^] 
studied particular solutions of the BD field equations with a charged perfect fluid source in Dicke's conformally 
transformed units, van der Bergh analysed the case of a ST spacetime containing an electrostatic field, with 
traceless energy-momentum tensor, and found both naked time-like singularities and non-singular solutions. Reddy 
and Reddy & Rao also obtained spherically symmetric, static, conformally flat solutions of the BD electro- 
vacuum field equations. Neither of these authors, however, have considered the case of massive particles bounded to 
form a star, and the existing literature does not provide an analysis of the charged boson star system as we do here. 

While charged stars, either in GR or in a ST theory, are interesting objects from a theoretical standpoint, until last 
year it was thought unlikely that they would play a significant astrophysical role. The argument to support this view 
was that of selective accretion: any charged object will naturally accrete matter of opposite sign and ultimately become 
neutral, as would, for example, a Kerr-Newmann black hole. However, a recent study by Punsly established that 
a charged rotating black hole will be able not only to preserve its charge, but will also generate observable gamma 
rays. In an extremely simplified view of his model, he considered the black hole to be surrounded by a ring of 
opposite charge, but otherwise isolated in interstellar space. The ring, while neutralizing the whole system as seen 
by an observer at spatial infinity, would allow the black hole to keep its charge. Magnetic flux lines would enter the 
hole through its axis of symmetry and pair creation would produce gamma rays. This process may be important in 
understanding the nature of some of the unidentified galactic gamma ray sources recently discovered by the NASA 
Compton satellite, in its EGRET experiment psj . 

The already initiated program of searching for observational signals of boson stars has lead us to some significant 
analogies with black holes |2g]. Therefore, it is natural to expect that research into the subject of charged boson 
stars would extend this analogy. For instance, it may be possible that a rotating charged boson star may replace the 
central engine of Punsly's mechanism, although a careful analysis should be made to prove or disprove this statement 
(something which is far from the purpose of the present work). It is worth noting, however, that the more realistic 
models of (uncharged) rotating boson stars have been investigated and solutions were shown to exist ||2^. There is 
no reason to suppose that rotating charged boson star solutions do not also exist. All in all, it is encouraging that 
charged objects are found to produce testable effects from which their astrophysical signatures may be extracted: the 
presence of charge allows a wider range of phenomena to be observable. In the context of our present work, boson 
stars, either charged or uncharged, are relativistic constructs that may help us understand the influence of a different 
gravitational theory and its cosmology on the astrophysical world. 

In a ST theory, at the level of the action, Newton's gravitational constant G is replaced by a field G* = G/i/). 
In a cosmological setting G* will evolve with time, giving rise to the twin phenomena of gravitational memory and 
gravitational evolution p^ . Consider a compact object in equilibrium that forms at some time ti at which the 
gravitational coupling strength is G*. The value of will determine to some extent the structure of the object. At 
some later cosmological time, G* will have evolved to some new value G2 and the object may either "remember" 
its local value of G* or evolve to some new configuration in which its local value of G* matches the cosmological 
value Gj. In actual fact, one would expect the behaviour of the object to lie somewhere between these two extremes. 
Boson stars seem to be ideal candidates for the study of these phenomena: they are relatively simple objects whose 
equilibrium solutions are easy to find, they are non-singular and they are fully relativistic. Gravitational memory and 
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gravitational evolution in boson stars were first analysed by Torres, Liddle & Schunck ||2^ , who constructed sequences 
of static configurations at different cosmological times. One can also consider the evolution of other compact objects, 
such as the white dwarfs studied by Garci'a-Berro et al. |Q and Benvenuto et al. , who found that the variation 
of G* produced observable changes in the stars' luminosities. 

The bosons that make up a boson star are assumed to have identical mass m and, as we shall show in detail in the 
following Section, m and the Planck mass Mpi serve as a scale against which all physical properties of the star may 
be measured. Here we shall briefly outline the range of physical magnitudes possible for given choices of m. For a 
mini-soliton star (made of uncharged non self- interacting bosons), the choice m = 30 GeV implies that the mass of the 
star is of order 10^° kg, which is 20 orders of magnitude smaller that the solar mass Mq, while its radius is of order 
IQ-i^ m. This gives a density of around 10^® times that of a neutron star Q and these objects are obviously highly 
relativistic. Only when the boson mass is reduced to around 10~^°eV does the star resemble a more conventional 
stellar object. As we have mentioned above, including a self-interaction potential V(ip) can substantially alter these 
figures. The importance of the self-interaction term is measured by the ratio V{ip)/{m'^\^\'^)-, which for the quartic 
self- interaction considered by Colpi et al. [D is approximately XM'^i/w?, where A is a constant and one assumes that 
IV'I ~ Mpi. In this case, the self-interaction can only be neglected if A is exceedingly small: A <C m^/A/^j = 6.7 x 10^"^^ 
GeV~^m^. Rewriting the potential as the dimcnsionlcss number A XM'^J {^irra^), the equilibrium solutions can 
be parameterised by A and for A 1, the boson star masses scale as A^/^M^j/m ~ A^/^M^j/m^ ~ A^/^Afch, where 
Mch is the Chandrasekhar mass (approximately IMq). For instance, if m = 1 GeV, the star has a mass of 1.6\/AMq 
and a radius of 6\/Akni. If, instead, we take m = 1 MeV, then the mass becomes lO^-x/AM© and the radius becomes 
10^-\/A km, which is similar to the Sun's radius but encompasses one million solar masses. For non-topological soliton 
stars, the mass was found to scale as M^^/rn? and one can construct soliton stars of very large mass and radius 
(the latter quantity being of the order of a few light years). The fact that both m and A are free parameters allows 
one to construct boson star solutions with physical properties whose values range over all astrophysically interesting 
magnitudes. With the inclusion of the gauge charge, the range of physical magnitudes will depend upon a third 



independent parameter. We shall discuss this in more detail when we analyse numerical solutions in Section VI. 
The rest of this work is organized as follows. Section |l| introduces the basic formalism, both for gravity and the 



stellar model. An analysis of the charge, the mass definition and the asymptotic limit is made in Section III. In 



Section |V| we discuss the boundary values for all fields and provide scaling relations valid for our field equations, 
which we shall later use. We then analytically prove, in Section ^ that there exist a maximum charge for weak 
field configurations (characterised by small central densities) which differs from the GR result (we shall also explicitly 
derive this last result as the correct limit of the ST case). In Section ^ we numerically analyse strong field solutions 
for both BD gravity and a representative power law ST theory. These latter solutions demonstrate how ST boson 



stars may produce observable ST effects at the current epoch. Section VII is devoted to the discussion of gravitational 



memory and evolution of boson stars, both charged and uncharged. We give our final comments in Section VIII . 



II. FORMALISM 

The system we study is that of a complex scalar field i/" with a C/(l) charge in the framework of ST gravity. Its 
action is given by 



IQtt 



(1) 



where g — Det^g^^i,), R is the scalar curvature, u is the coupling function and is the matter Lagrangian. We take 
Lm to be 

Lm = -lg''-(D^'^){D,yj) - im^?- ]k^W - 7^^^"^^^- (2) 
z z 4 ion 

where 

Ffiij dpijAiy djjj\^^ (3) 
D^,=d^+ieA^, (4) 
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Af_i is the gauge field and the over-bar denotes complex conjugation. This Lagrangian is invariant under a gauge 
transformation of the U{1) group. This implies the existence of a conserved current 

J'' = \/^5^" [«e - i'd.i^) - 2e^A,H^ , d^J'' = 0. (5) 

Varying the action with respect to g'^'^ and we obtain the field equations 



1 



and 



where 
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2w + 3 



SttT- 



duj 



.■,v 9fj.i/ 



An 



F '^F 



1 

8^ 



is the energy-momentum tensor for the matter fields and T is its trace. 

Varying the action with respect to the mater fields V', 4' ^'Hd we obtain the boson wave equations 

g^^VaDpi; = (to^ + A^i/')V' - ieg^'^AaD^i;, g'^'^Vj^fi'^ = (m^ + A^?A)V^ + ieg'^'^ Aj3 p'^ 

and the Maxwell equation 

where 

1 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



We consider static, spherically symmetric solutions and write the line element in the Schwarzschild form 

ds^ = -B{r)dt^ + A{r)dr^ + r^d^^, (12) 

where dVi^ is the line element of the unit 2-sphere. We look for solutions of minimum energy. It has been shown 0| 
that for spherically symmetric systems of uncharged bosons, since Lm is universally coupled to the curvature, the 
form of ip compatible with this requirement is 



•^{r.t) = x('')exp [iwt], 



(13) 



where 137 is a real positive constant. To prove this, one minimises the total energy of the boson star subject to the 
constraint that the particle number (which we shall define below) be conserved. The parameter ■uj plays the role of 
a Lagrange multiplier and this result is independent of the form of the boson self-interaction V{<j)). This result also 
holds when one includes a U{1) charge and so we take Eq. |l^to be valid in the present work. 

Equation (|l3| ) is also consistent with the assumption that is static. The scalar field inherits the symmetries 
of the line clement and is a function of r, while the gauge field may be chosen such that we have only electric charges 
and no magnetic ones, so we set 



A^, = (A(r-), 0,0,0). 



(14) 



Before we explicitly write down the full set of equations of structure of the star, we shall introduce some dimensionless 
quantities which will simplify the comparison with the previous work cited in Section ||. Using the boson mass m and 
the present value of the Planck mass Mpi :— yl/G to define length and energy scales, the new quantities are: 



n{r) 



m + cAq 



(15) 
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and 



for the boson sector, 



and 



a = \/47r' 



M 



pi 



pi 



eMpi 



(16) 



(17) 



(18) 



(19) 



for the gauge field, and 



for the Brans-Dicke field. We also introduce a rescaled radial coordinate x defined by 

X := mr. 

In terms of these quantities, the independent components of the field equations (^, (0), ^ and ( |lO|) are 
• The generalized Einstein equations: 

A' = Ax -rrz :rr- + —{ (2cj + 5) + 2cj - 1 ' ' ^ ' 



$(2w + 3) $(2w + 3) V B 



2B$ 2$2 $(2w + 3) Vrf*/ S$ a; 



2$(2w + 3) 
duj\ C"2 1 ^ 



(20) 



(21) 



(22) 



2$ + a;$' 



2Ba'2 + 2a'A{n^ - B) - rx^ABA - — + —$'2 - 2C'2 + ^(A - 1) 

a; $ a; 



(23) 



• The scalar field equation: 



\2A 2B X 



2A 
2uj + 3 



B J A 



1 duj 
2w + 3d¥' 



/2 



(24) 



The boson field equation: 



^ [2A-2B-x^-^ 



(25) 



• The Maxwell equation: 



(26) 



Here a prime denotes d/dx. We mainly consider two forms of lu: the Brans Dicke coupling in which uj is constant and 
a power law coupling, where 2w + 3 = 4n/3$". 

Equations ( p2[ ) to ( |2^ ) reduce to the general relativistic ones Q when uj — > oo, <i> ^ 1 and dio/d^ —i- 0. Note, 
however, that this is not a simple issue whenever the trace of the matter energy- momentum tensor vanishes [ p3[ . This 
is not the case here, although we still verify the general relativistic limit for large values of w. 
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III. CHARGE, MASS AND THE ASYMPTOTIC LIMIT 



We shall first look for solutions that are asymptotically flat in the sense that, in the limit x 



A' 



l + 0[ — 



Oil- 



Oil- 



o — 



(27) 



where n > 1 and the subscript 'cx)' denotes values at space-like infinity. We also require a to vanish at least as fast 
as so that the boson matter is localised and the solutions describe bound eigenstates of ip. In fact, as we shall 
show below, the conditions above imply that cr falls off exponentially as a; ^ oo. In addition, the field equations ( |2^ ) 
to (|2^) are invariant under the rescaling C^C + c, w—>w — qc where c is a constant, and we use this remaining 
gauge freedom to set 



0. 



(28) 



This implies that Q^o — ^/m. 

Given the asymptotic flatness conditions, we may deflne expressions for the total charge and mass (or energy) of 
the boson stars. From the conserved current equation (^ one may derive a conserved (time independent) charge 



A 

fl\ — cr'^x^ dx. 
V B 



(29) 



This quantity may be interpreted as the star's total electrical charge per unit boson mass. Since all of the bosons are 
assumed to carry identical charge per unit mass q, we may also define a conserved particle number 



TV 



Q\ — cr^x"^ dx. 



(30) 



Due to the fact that we are using the boson mass m as the mass scale, the numerical value of N is also equal to the 
value of the 'rest mass' of the star, a term we use to describe the sum of the masses of the bosons measured in some 
non-gravitational way. Then the particle number is measured in units of M^j/rn^ while the rest mass is measured in 
units of Mpi/m. We also define a Newtonian mass 



N 



N 



(31) 



which measures the active gravitational mass of the star whose bosons are dispersed to infinity. 
Integrating the Maxwell equation (pq) we find an alternative expression for the total charge: 



d-x 



dC 



fAB dx 



dx = lim ( x'^ —— 

dx 



(32) 



where we have use the asymptotic flatness conditions to derive the second equation. In an analogous way, we deflne 
the 'scalar charge' 



5* 



lim 

X — ^oo 



dx 



(33) 



Note that this is not in general a conserved quantity as there is no internal symmetry associated with $ in the action 
(|l|). We introduce it here since it is included in the expression for the mass of the star which we shall give below. 
Integrating the wave equation (E4h gives an alternative expression for the scalar charge: 



S = 







2A 2B 



2A 



B 



1 doj 
2^7+3 



2u} + i 

The generalisation of the Schwarzschild mass is deflned implicitly by the familiar relation 

-1 



dx. 



A:= l- 



2M{x) 



(34) 



(35) 
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Note that this equation does not exphcitly include the scalar field $: here the gravitational coupling strength G* 
G/^ has been factored into M. Using Eq. the Schwarzschild mass may be written as the integral 



M(x) 



J2 

dx- 



cLu 



+ 



(2w + 5) + 2w - 1 



A$(2w + 3) 2$(2w + 3) 2ylS$ 



(36) 



One can show that Moo '■— limx^oo A/ = Madm, where Madm is the ADM mass of the star. Note that the 
integrand in Eq. (|3^) contains non-positive definite terms, so that M' and M themselves may be negative for weak 
field solutions. This is because the scalar field terms on the right hand side of Eq. (||) may cause G^^ to violate the 
dominant energy condition. This fact has been exploited in the construction of ST wormhole solutions for which the 
matter energy-momentum tensor satisfies all of the energy conditions and even when there is no matter tensor at 
all [H- ^ 

As shown by Lee in a usually uncited paper ||3^ , the correct definition of mass for an isolated source in BD gravity 
is the Tensor mass 



M 



S 



ADM 



2$, 



(37) 



This definition is also appropriate for more general ST theories. From the Tensor and Newtonian masses, we define 
the binding energy 



£ := Mt - Mat 



(38) 



and the fractional binding energy 



Mt - M^ 



N 



Mn 



(39) 



which measures the binding energy per unit boson mass per boson. A necessary condition for dynamical stability is 
that f < 0, which implies that B <Q. 

There are several reasons for preferring Mt over Madm as the correct definition of the mass of the star in ST 
gravity. Lee |3^ has shown that for a localised source that emits gravitational radiation, the time rate of change 
of Mt evaluated at future null infinity is non-positive definite (in other words, gravitational wave emission can only 
reduce the Tensor mass of the source). This is not true for Madm- Creighton and Mann have given a quasi- local 
definition of mass in ST gravity based on a rigorous consideration of the boundary terms both in the variation of 
the Lagrangian and in the formulation of the Hamiltonian. It turns out in either case that the expression for the 
mass is given by the usual integral of the extrinsic curvature on the boundary of the region under consideration, plus 
an additional term analogous to the scalar charge defined above. For an asymptotically fiat spacetime, and in the 
appropriate limit, one can show that their mass tends to the Tensor mass given above. In addition, one can show 
that, at least for uncharged boson stars in BD theory, aiw solution with extremal Tensor mass is also one of extremal 
Newtonian mass (and hence extremal particle number) [|38| . This result is easily extended to more general ST theories 
and is important when one uses catastrophe theory to determine the stability of a boson star |39| : a plot of mass 
against particle number for a set of solutions should show "cusps" at the extremal values of the particle number. The 
first of these cusps marks the point at which the solutions become dynamically unstable. 

Expanding the metric and matter field equations in powers of 1 /x about x — oo one can write the field equations 
in a linearised form, similar to that given in Section ^ for the weak field limit. Solving these equations and using the 
boundary conditions (p^) and (Eq), we find that the scalar field and vector potential have the asymptotic form 



Q / 1 

X \x^ 



X V 



(40) 



while the line element in this limit is 



ds^ = - 



1 - 



2Mk 



x^ 



dt^ 



'2.MADM I g ( 1 

X \x'^ 



dx' 



(41) 



Here 
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Mk Madm — -T — 



(42) 



is the Keplerian mass of the star. Note that Eq. ( [4Cl| ) can be obtained directly from Eqs. ( p2D and (|33 

Substituting the metric components appearing in the asymptotic form of the line element into the boson wave 
equation (|5|) gives 



X 

This equation has the solution 



a 1 



2Madm \ 
X ) 



1 



2M 



K 



- 1 



x"^ 



where the constants k and h are given by 



n = ^i-ni^, h = -l + - 

K 



i + o{- 

x 



{2nl - i)Mt 



2$. 



(43) 



(44) 



(45) 



Hence the boson field falls off exponentially with x in the asymptotic region. The quantity k may be interpreted as 
the reciprocal radius of the star and must be real. This implies that fioo < 1. 

We conclude this section with a note on the physical interpretation of Mk and Mt and the choice of the physical 
frame. The Keplerian mass is the active gravitational mass measured by a non self-gravitating test particle in a 
circular orbit at space-like infinity about the star, while the Tensor mass is the active gravitational mass measured by 
a 'test black hole' in a similar orbit |40j. By test black hole we mean an object whose mass is negligible compared to 
the mass of the boson star and is made up entirely of gravitational binding energy. Its orbit is geodesic in the Einstein 
frame, whose metric g^i, is related to the physical metric by the conformal transformation 5^1, — ^g^^. Calculating 
the angular velocity of a circular orbit in the Einstein frame, using Kepler's third law and transforming the result 
back into the Jordan frame, we find that My is given by 



Mt = lim 



r2($B' + B$') 



2$ + rcj)' 



(46) 



The expression for Mt in the Einstein frame shows that this quantity is non-negative. 

The ADM mass in the Jordan frame cannot be interpreted in terms of particle orbits. In the limit of large cj, $ 
is only weakly coupled to the curvature and S is small compared with Madm- Hence, in this limit the differences 
between Madm, Mk and Mt are negligible and in the exact GR theory all three definitions are the same. 

Recently, Faraoni, Gunzig & Nardone have argued that the Einstein frame should be considered to be the 
physical frame, primarily because the Jordan frame ADM mass is non-positive definite, even when the matter energy- 
momentum tensor T^^, satisfies all of the energy conditions. However, we believe that this reasoning is based upon the 
wrong choice of mass definition: the physical mass in the Jordan frame is Mt and one can show that Mt is positive 
in the following way. After a conformal transformation into the Einstein frame, one can show that the new Einstein 
tensor will obey the dominant energy condition if the Jordan frame matter energy-momentum tensor T^^ obeys the 
dominant energy condition. Then the Einstein frame ADM mass will be non-negative, by the usual positive energy 
theorems. However, by analysing the effect of a conformal transformation on the definition of Madm, one can show 
that the Einstein frame ADM mass is identical to the Jordan frame Tensor mass, up to a (positive) factor of $00 
p2[ . Recently Santiago & Silbergleit |Q have rewritten the Jordan frame field equations of ST gravity by defining 
a new connection that isolates the dynamical degrees of freedom of the metric gravitational field. Their formalism is 
algebraically equivalent to making a conformal transformation into the Einstein frame, so it not surprising that their 
result corroborates with the proof outlined above. The scalar field contribution to the total energy-momentum of the 
spacetime can now be made non-negative, contrary to the previous claim of Faraoni et al. Here we simply use the 
conformal relationship between the two frames as a computational aid, and give our results in the physical Jordan 
frame. 



IV. BOUNDARY VALUES AND SCALING RELATIONS 



In addition to requiring asymptotic flatness, we impose the boundary conditions 
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^0 = 1, a'^^^'^ = C'^ = 0, (47) 

where the subscript '0' denotes values at a; = 0. These conditions ensure that the solutions are regular at the origin 
and, combined with the asymptotic flatness conditions (|2^), they make the field equations eigenvalue equations for 
Tu which automatically lead to the asymptotic conditions cr^ = B'^ = = = 0, A^o = 1. The values of $00 
and Boo are not determined by the field equations, and we next show how they may each be rescaled to take on any 
value. 



The field equations (22) to (^6|) are invariant under the global rescaling 

B^-fB, vo^^w, C^^C, n^^n, (48) 
where 7 is some constant. Equation (Ha) simply rescales the magnitude of the time-like Killing vector field, and 



leaves all solutions physically unchanged. With the gauge freedom allowed by Eq. (48), the field equations become 
eigenvalue equations for Bq with il arbitrary and their solution in general leads to a value of B^o 7^ 1. We use Eq. 
( ^ ) to set Boo — 1 for each solution. 

The equations of motion are also invariant under the global rescaling 

a^ka, C-^kC, q^r, A (49) 

k fc^ 

where k is some constant, provided w(<i>) is held invariant. For a general ST theory, this will require a change in the 
functional form of uj. 

Equation (|4^) leaves the ratio S/ ^00 and the ADM, Keplerian, Tensor and Newtonian masses invariant and rescales 
the total charge and particle number as 

kQ, N ^ k^N. (50) 

Physically, the rescaling exchanges gravitating matter (measured by the number of bosons) for gravitational field 
energy (which depends qualitatively upon the gravitational coupling strength G* = G/^) in such a way as to leave 
the total mass of the boson star invariant. 
Consider now the set of solutions 

5(cto;$oo,'?, A) = {Af(CTo;$oo,g,A), iV(o-o;$oo,g,A), Q(cto;*oo,9,A)} (51) 



parameterised by ctq for fixed $00, q and A, where M represents all of the mass values discussed in Section IV. Under 
the rescaling $00 k^^oc the mass curves are re-parameterised as 

M{ao; $00, g, A) -> M{kao; e<^oo,k-\, k'^A) (52) 

while the particle number and charge curves are both re-parameterised and rescaled as 

N{ao; $00, 9, A) -> k^N{kao; k^<i>o., k'^q, k'^A) (53) 

<9(c^o; $00, g, A) ^ kQ{kao; k'^<i>oo, k^^q, k^^A). (54) 
Hence, by applying the rescaling (^9|), S may be used to generate a new physically distinct set of solutions 

5(CTo;fc^^oo,g,A) = {M(cro;^^$oo,g, A), iV(cro;fc^«'oo,g,A), g(cro;fc^«'oo,g,A)} 

= {M(fc- Vo; ^oo.kq, fc^A), k'Nik-^ao; ^00, kq, k^A), A:Q(fc-Vo; $00, A^g, fc^A)}. (55) 

Thus, Eq. ( ^ ) may be used to compare boson stars in cosmological settings with different values of <i>oo. 

The existence of this scaling relation means that we need to vary at most two of the parameters $00, A and q in 
order to investigate completely the properties of a solution in any given theory. For example, by fixing $00 to some 
constant, the complete set of solutions may be generated by the subset for which only A and q are varied. In addition, 
if either of these two quantities is zero, we need only vary the other to build up a complete set of solutions. This will 
be of use when we discuss our numerical results in Section There we shall be primarily concerned with stars for 
which A = and we shall show complete sets of solutions for several theories by only considering different values of q. 

Finally, we comment on another possible choice of mass and how it is affected by Eq. (^|). Torres, Schunck & 
Liddle define the Schwarzschild mass by the relation 
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Note that A{x) is the same metric potential both here and in Ref. . Comparing this expression with Eq. ( p5[ ) we 
have M* — ^ooM, which implies that M* rescales by a factor of under Eq. (Eoh with q ~ 0. This will be important 



when we analyse the behaviour of AI* in a star undergoing gravitational evolution (see Section VII). They then define 



the binding energy by using = limaj^oo A/* together with the rest mass (as we have defined it) to give 

£* = A/;, - N, (57) 

which also rescales by a factor of under Eq. (^9|). Their choice is such that the binding energy of any solution has 
the same sign when compared with ours and the fractional binding energy (defined a,s B — £* /N in p^ ) is identical 
for both papers. Thus, if 



then 



m;, - iV < (58) 



M* M* N 

M^-Mm = -r^-Miv = ^-:;-<0 (59) 



and 

N ~ Mn ' 



(60) 



Here we have used the ADM mass (as we have defined it) instead of the Tensor mass since Torres et al. use the former 
in [Q. As they point out, for the large values of uj they consider, the numerical difference between Madm and Mt 
is negligible. 

Since M* and the masses we define in this paper are measured in the same units, it is apparent that they measure 
physically different quantities. As mentioned in the text after Eq. (|35|), the masses we define implicitly include a 
factor of l/^oo- This means that Mt (or Madm when uj is very large) measures the total gravitational energy and 
Mk measures the active gravitational mass as seen by an orbiting test particle with negligible self-gravity. Our mass 
definitions correspond to the product GAA of Newtonian theory, where is the Newtonian rest mass. The Newtonian 
analogue of Eq. (^) is ~> k^M, G — > G/fc^ which leaves GM invariant: the rescaling swaps energy between the 
matter and the gravitational field. In the weak field limit and when w is large, describes the rest mass of the 
star in which the gravitational coupling strength l/$oo has been factored out (see the discussion after Eq. ( |83| ) in 
the following Section). Thus M^^ roughly corresponds to M.. However, for strong field solutions one cannot easily 
interpret , since one cannot decouple the gravitational field strength from the total energy. 



V. WEAK FIELD LIMIT AND THE LIMITING CHARGE 



In this Section we discuss the weak field limit of the solutions which we define as the limit in which co is small but 
non-zero. The analysis that follows is in part a generalisation of the approach followed by Kaup |Q when dealing 
with uncharged boson stars in GR and involves expanding the solutions in power series about fiat spacetime. We 
start by defining a function n(x) which has the boundary value Hq = 1 and the functional form given implicitly by 

a = eH. (61) 

The constant parameter e :— (To/Ho measures the degree to which the solution differs from flat spacetime and we shall 
use it as the expansion parameter in the analysis. We also define a rescaled radial coordinate by a :— y/sx. Hence, 
H is now considered to be a function of a. Note that, by definition, H(a) has the same asymptotic form as a and 
decreases exponentially as a — > oo. 

The field variables and energy parameter we expand in powers of e about flat spacetime. We implicitly define new 
field variables a, f3, S and F by the first order expansions 

A = l + sa{a) + Oie'^), (62) 



10 



B = l + e0{a) + O{s^), 



(63) 



and 



$ = $oo+eS(a) + 0(£2), (64) 



C^tr(a) + 0{e') (65) 



n = l-^-E+^-qT{a) + 0{s^), (66) 



where $oo is the value of the scalar field at space-like infinity and i? is a constant for each solution. We assume that 
the scalar field coupling parameter has the form 

u;^Lu^ + eW{a) + Oie^), ^^A{a)+0{e), (67) 

where Woo is a constant, A = lime^o (^) ^-nd we have written the second term in the expansion for uj explicitly in 
terms of a. These assumptions restrict the applicability of the following analysis to a subset of all possible scalar-tensor 
theories. However, Eqs. ( |67|) are satisfied for both BD theory and the power law theory. In the latter case we have 

-=^('J'oo+eE + 0(s^))"-| (68) 

Expanding this expression binomially, we see that the power law theory has W — 2n^E/(3<i>oo) and A = 2n^/(3<l>oo)- 
In terms of the new variables the field equations (^2|) to ( p6| ) become 

4n2a2(woo + l) , . 
^""^ ^ a>^(2..^ + 3) +^(^)' ^''^ 

P' = -~^f + Oie), (70) 

$oo(2Woo+3) 



n')' = n(i? + /3-gr) + o(£), (72) 



and 

(a^T')' ^2qll^a^ + 0{e), (74) 



where in this Section a prime denotes d/da. We have derived Eq. (|7l| ) by differentiating Eq. (|70| ) and substituting in 
Eqs. (|69| ) and (fz^). Equations ( ^2|) to (|65| ) taken together with the asymptotic conditions (|27|) and (^ ) imply that 
we must impose the boundary conditions 

oioo = /3oo = Too = Soo = 0. (75) 

Equations (|7l|), (^3|) and (^4|) imply that both (3 and F are strictly increasing functions of a while S is strictly 
decreasing. Combining these properties with the boundary conditions above we have /? < 0, F < and S > 0, where 
the equalities are only true in the limit a ~* oo. From Eq. (|6^), the requirement that iloo — nj/m < 1 and the 
vanishing of Fqo imply that -B > 0, where the equality is only true when 11 = 0. 
Combining Eqs. (||) and (^) we have 
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Q^\^e lim ia^^ = ^fe H qU^a^ da + O {e'^') := V^qAf + O {e'^') . 
a^oo \^ daj Jq 



(76) 



where we have defined 



Af: 



a^n^ da. 



(77) 



The quantity ^/eM is the rest mass (or particle number) to lowest order in e. From the definition (^) the Newtonian 
mass in the weak field limit is given by 



(78) 



From Eqs. (|33|) and ([73^, the scalar charge is given by 



S — \fe lim a 



da 1 



2uj. 



2UJry 



3 



(79) 



where we have used Eq. (|7^) to obtain the final equality. Integrating Eq. (^9|) and using the definition (^5|), we obtain 
the weak field ADM mass 



$oo(2t^oo +3) 



(80) 



where again we have used Eq. (|7^). Combining Eqs. ( [79[ ) and ( |80| ) with the definitions ( pTf ) and (42) gives the weak 
field tensor mass 



(81) 



and the weak field Keplerian mass 



$oo(2l^oo+3) ^ ^ 



(82) 



Equations (^Tj) and (|2|) show that the Keplerian and Tensor masses may be written as the products Mk = ^/sGk-^ 
and Mt — y/eGx.^, where the coupling strengths Gk and Gt are given by 



Gk 



$oo(2cJoo+3)' 



Gt 



(83) 



Hence, the weak field solutions are Newtonian in the sense that the Keplerian and Tensor masses, which are the 
active masses to which test particles respond, may be decomposed as products of the star's rest mass and a coupling 
strength. In addition we have Mt = Mn to lowest order in e, which motivates the labelling of the quantity defined 
in Eqn. ( |3l| ) as the Newtonian mass. 

To calculate the fractional binding energy, we need to write Mt and Mn up to order e^^^. Rewriting the integrals 
(53) and (Isl) in terms of the weak field variables, one can show that the ADM mass and scalar charge are given by 



M 



ADM 



$oo(2Woo +3) 



a 



1) 



F'^ 



2(2tJoo + 3) 8 



4$. 



2(2cjoo + 3) (2cjoo + 3) V 2t^oc + 3 



2W 21] 

+ (2cjoo + 5){qT -E-P)- — {uJco + l) 



and 



'2^U 
2tJoo + 3 



2n'- 



Ay/2 yl 



+ 



2n^ 

(2c^oo + 3) 



qT-E- f3~a 



2cjoo + 3 2cjoo + 3 2 
2T4^ 



(2c<Joo + 3) 



da + O (£=/^) 



da + 0(£^^^) 



(84) 



(85) 
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Similarly, rewriting Eq. (p9|) and using the definition (BG) we have 



Mn = 



-3/2 



2$c 



Combining Eqs. (pTl), (BJ), (|85|) and (Bq) with the expression (|39[) for the fractional binding energy we have 



^ Jo 



E'^w^ EV n'2 r'2 ^-a(2tJoo + 1) 4S(tJoo + l) 



4$, 



2wo 



$oo(2Woo +3) 



(86) 



(87) 



We need to simplify this expression considerably. To do this we integrate Eq. ( p7[ ) by parts in order to write each term 
as some numerical factor times the same positive definite integral. We carry out each partial integration separately 
below, noting that all surface terms vanish due to the asymptotic conditions (|7^) and the fact that all of the fields 
are finite at a = 0. 

Using the wave equation ( [t^ ) one can show that 



j*oo poo 

/ a^n'^da^- a^n^{E + /3~qr)da. 
Jo Jo 



Integrating the product aW{a^Il'y by parts and using Eq. (|8q ) one can also show that 



/•CXD /"OO 

/ a^U'^da^^ a^U^iP' -qV) da. 
Jo Jo 



Integrating a^S'a' by parts and using Eqs. (Foh and ([73|) we have 



a^Ti'a' da 



°° 2n2a3 



2E' 



da. 



(89) 



(90) 



Integrating both a^E'^ and the product aE'(a^E')' by parts, using the wave equation ( |7^ ) and combining the results 
one can show that 



a^E'^ da = - 



4a3n2E' 



da = 



da. 



Similarly, one can also show that 



/•CxD nOO POO 

/ a^r'^ da= Aa^lfqT' da = - 2a'^n'^qr da. 
Jo Jo Jo 



Using Eq. (B3) we have 



a^li^a da 



CO / OV' 

3tt2 a' , 



$0 



da. 



Substituting Eqs. (p9|) to (M) into our expression for the fractional binding energy mm gives 



or' 1 
- — \ , 

4 2woo + 3 V * 



:^(^oo + 2)-^(2c^^ + l) 



da + O (e^) 



To simplify this expression further, we use Eq. ( fz^ to obtain 

/•OC /"OO /"OO 

/ a^tfr' = / atf (a^r') da = / 2qXan2 da, 
Jo Jo Jo 



where we have defined a new function 



X{a) := / d^jf dd. 
Jo 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 



Similarly, from Eqs. ([Tl]) and ([73) we have 
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a^tf^' da = / ^ rfa 97 





and 

a^tfE' da = - / \ da. (98) 

Jo 2ajoo + 3 

The quantity iP is positive or zero over the interval of the integration which implies that X is an increasing positive 
function of a that vanishes at a = 0. The exponential decrease of 11 in the limit a — > oo implies that X has a finite 



limit Xoo- Hence the integrals (95), (|97| ) and ( |98|) are finite and substituting them into Eq. (|9J) gives 



J\f\2 $oo(2woo + 3) 



/ aXn^ da + O (e^) . (99) 
Jo 



This result is independent of A and is true for any form of the coupling parameter a;($) that satisfies Eq. (|C 
Equation (p9|) implies that the fractional binding energy is negative provided 



,<,^.A^M-^^§^y (100) 

For the star to be stable B must be negative. Hence Eq. ( |lOCl| ) places an upper limit on the charge to mass ratio of the 
bosons that form a stable star in the weak field limit. For a star with q > qmax the Coulomb repulsion dominates over 
the gravitational attraction and no stable solution exists. Note that, as one would intuitively expect, Qmax oc \/l/$oo 
for fixed ui, so that the greater the gravitational coupling strength, the greater the amount of charge that can be bound 



into a stable object. In Section VI we shall numerically examine how Qmax is modified for strong field solutions. 

For finite lo we have Qmax > 1 for $oo = 1 which implies that in ST gravity it is possible to construct a stable 
star from bosons whose charge is greater than that allowed by GR. In the GR limit, where a; — > oo, Eq. ( |100| ) gives 
Qmax = 1, which is the upper limit on q given by Jetzer and van der Bij by just comparing the electrostatic and 
electromagnetic forces (note that our unit of charge differs from theirs by a factor of a/2)- Hence, the above 
derivation serves also as a detailed proof of the maximum charge in General Relativity. 



VI. NUMERICAL SOLUTIONS 

In this Section we discuss strong field solutions of the field equations (^2|) to ( p6|) for both BD theory and the power 
law ST theory. In their full form, the equations are complicated enough to require numerical integration. To generate 
each solution in the set 5, the numerical routine that finds the eigenvalue must be nested inside a second iterative 
routine that searches for the the value of $o that gives a value of <I>oo common to all solutions in S. When generating 
each solution, the numerical routine halts at the value of x at which a becomes too small for the code to proceed 
further, typically at a value a ~ lO"^*'. Then the Tensor mass is calculated using Eq (^) with the limit replaced by 
values of x near to the termination point of the routine and a power series approximation about a; = oo is used to 
estimate Mt for the solution. A similar power series method is used to calculate N and iS, based on the behaviour of 
$(2:) and the integral of dN/dx out to the termination point. 

For the sake of ease of comparison, we have chosen $00 — 1 for all of the solution sets we discuss here, and this 
choice implies that Mjv and N are numerically identical for each solution. The sets of equilibrium solutions are shown 
in Figures 1 to 5. In each of these Figures we show Tensor mass and Newtonian mass curves as a function of do for 
several sets of solutions within the same theory, choosing a different value of q for each pair of curves. The results of 
the numerical integrations arc then used to construct Figures 6 to 9. We shall discuss each figure separately below. 

Figure 1 shows the mass curves for several sets of solutions in the oj — 500, A = BD theory, for which the weak field 
charge limit is qmax = 1-001. As one would expect from such a large value of lj, the curves in Figure 1 are virtually 
indistinguishable from the mass and particle number curves calculated for GR boson stars in [Q. As g approaches 
Qmax from below, the locations of the maxima in the Mt and Mjv curves shift to lower values of cto and the binding 
energies £ decrease in magnitude, while the gradients dAIx/dao and dM^/daQ of the curves near — Q increase. 
In addition, as q increases, the mass curves shift to successively higher values, since a star consisting of bosons with 
large q must generate a strong gravitational field to overcome its own Coulomb repulsion. This is achieved by an 
increase in N (hence an increase in Mjv) which leads to an increase in Mt- Since u is large for these solutions, $ 
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remains approximately homogeneous throughout each star in each solution set shown in the Figure. For example, in 
the maximum mass q = 0.99 solution, $0 and $00 differ by less that 1%. Hence the solutions in the Figure show no 
strong scalar-tensor gravitational effects. 

In Figure 2 we show mass curves for several sets of solutions in the ut ^ —1, A — BD theory. As mentioned in 
Section |, this choice of the coupling constant may be relevant to the study of boson stars in the very early Universe, 
although in this case our value of = 1 is physically unrealistic: to approximate a boson star embedded in the early 
(possibly string theory dominated) Universe, we should choose a value $00 < 1- However, from the scaling relation 
( p9| ) , rescaling the asymptotic value of <& merely rescales the horizontal axis of Figure 2 and leaves the form of the 
curves invariant (although for $00 7^ 1, Mn and N are no longer numerically equal and q must be rescaled). Hence 
the figure may be used as a template to generate cj = — 1 boson stars solutions with any value of $00 • 

For the w = — 1 BD theory, the weak field charge limit is qmax = V^- As Figure 2 shows, this limit does not extend 
to the strong field solutions (in contrast with the GR charged boson stars for which the limit qmax — 1 holds for all 
values of ctq)- The mass curves of the solution set for which qmax > q > qc, where qc is some critical value of q, diverge 
at finite values of ao, and the value of ctq at which this happens decreases with increasing q. The mass curves of the 
solution set satisfying q < qc do not diverge, although for the values of q presented here the value of ao at which the 
mass curves are maximal increases with q. We have numerically found that qc ~ 0.90 for the cu = —1 solutions. This 
value is approximate since the equations are very hard to integrate when q is close to qc- 

Physically, we may understand the reason for the divergent behaviour of the masses as follows. For each solution, 
the ratio <I>o/3'oo > 1 and this quantity increases with ao since $ is strongly coupled to the curvature. Hence, within 
any particular solution set satisfying q > q^ the active gravitational mass of each boson (measured qualitatively by the 
product TO$^^) decreases as ao increases. At some finite ao-, the number of bosons with charge q needed to generate a 
gravitational field sufficient to overcome their own Coulomb repulsion diverges. The results shown in Figure 2 imply 
that the limiting charge qmax is a decreasing function of ctq that has the weak field value in the limit ctq — > and 
asymptotes towards qc as ao 00. Solutions with q > qmaxi^ro) do not exist, just as in the GR case. 

Figure 3 shows sets of solutions for = 1 BD gravity with A = 0. The weak field charge limit in this case is 
^1.2. The coupling between $ and the curvature is weaker than it is in the uj = —1 theory and the mass curves show 
behavior intermediate between those in Figures 1 and 2. For values of q up to q « 0.9, the location of the maxima in 
the Mt and Mj^ curves decreases with increasing q. As q increases further, the location of the maxima starts to shift 
towards higher values of ao until the mass curves diverge at q — qc ^ 0.976. This value is smaller than the co = —1 
critical charge and the numerical solutions suggest that qc is also a decreasing function of uj for fixed ao that only 
approaches the GR charge limit as a; — s- cx). 

Figure 4 shows mass curves for the power law ST theory, where 2uj{^) + 3 = |n<i>". We have set n = 4 in these 
solutions; this parameter choice is in agreement with current observational constraints on this kind of ST theories p6| . 
These are based on primordial nucleosynthesis calculations, which for this theory are stronger than the constraints 
imposed by solar system tests. (In fact, weak field test do not impose any bound on the exponent of the power law. 
This is a consequence of the form of the matter era cosmological solutions j4^, which are such that w 00 and 
u!~^duj/d^ when t — > 00 for all n.) As with the BD solutions, the mass curves diverge for values of q greater 
than the critical charge, which for this theory is qc ~ 0.986. The n — 4 power law boson stars could conceivably be 
observed today, and it is interesting that their mass curves are very different from their GR counterparts. 

We next outline how one converts from the dimensionless charge unit we use here to the SI unit of charge. Rewriting 



the weak field relationship (IOC) as an equation that balances gravitational and electrostatic forces on a boson in the 



asymptotic region of the star, one finds that the charge e of each boson, measured in coulombs, is given by 



where e is the permitivity of the boson star. For a star of particle number N, the total charge Q, measured in 
coulombs, is 



For example, if we take m = 30Gey « 5.3 x 10~^® kg, then e ~ (4.5 x 10~^^)q C. For a boson star with N = q = 1, 
one has Q « 1 C. This charge is small. However, such a star will have a very small radius (approximately 10^^^ m) 
so that its charge density is very large. 

In all of these solutions we have found that, for a given N and q, the magnitude of the scalar charge S increases 
while the mass Mt decreases, as the coupling strength l/uj is increased. This implies that for a given number of 
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bosons of charge the magnitude of the binding energy £ increases with the coupUng strength. Thus the bosons 
become more tightly bound even though the mass is smaller in the strong coupling case. One can understand this 
apparent discrepancy as follows. For all of these solutions, the scalar charge is negative and, comparing Eqs. ( |3^ ) 
and (Ji^), this implies that Mk > Mt, where the difference between these two quantities increases with the coupling 
strength. It is even possible for Mk to exceed Mn for particularly strong couplings (see Q for an example). One can 
show that the fractional binding energies B are reasonably small for all of these solutions (even for the rather extreme 
uj — —1 solutions, for which B never exceeds 0.2), so that to a first order of approximation the bosons that make up the 
atmosphere of a given star may be treated as a cloud of non self-gravitating test particles. Here, the term atmosphere 
rather loosely describes the asymptotic region of the star which encloses most of the mass, electrical charge and scalar 



charge. However, as mentioned at the end of Section III, non self-gravitating test particles respond to the Keplerian 
mass. Hence, in Newtonian terms, the gravitational force experienced by the bosons in the atmosphere tends to 
increase as the coupling strength increases, even though the energy of the star decreases. For a star with q < qc, the 
bosons in the central region of the star experience much the same conditions. The effects described here are entirely 
separate to the weakening of the gravitational field strength within the star caused by the increase of the ratio $o / 'i'oo 
with coupling strength. Although this reduces the Keplerian mass, in a star with q < qc this effect is swamped by the 
increase in the scalar charge. For stars with q > q^ the weakening of the gravitational field strength dominates over 
the increase in scalar charge and the divergent behaviour discussed above occurs. The fact that the scalar charge is 
large for strongly coupled stars will be important when we discuss the scalarization phenomenon below. 

Finally we consider the effect of introducing a quartic self-interaction term into the matter Lagrangian. For the 
GR solutions, the inclusion of the sclf-intcraction increases the mass of the stars but does not alter the value of qmax 
0. We have found that for ST theories, the inclusion of the A term not only increases the mass but also slightly 
decreases the charge limit of the strong field solutions. To illustrate this, Figure 5 shows mass curves for the w = 1 
BD theory in the limit A oo. The field equations in this limit are obtained by making the change of variables 



a ^ (J* = crA^^/^ and x ^ x* — xA^^'^, substituting these new variables into the field equations ( p^ ) to (26), and 
taking the limit A ^ oo (see and ||l^ for details). The resulting field equations are identical to Eqs. (22|) to 
( p6| ) except that the a'^ terms vanish, a'^A is replaced by cr*^, a and x are replaced by a* and x* respectively, and 
Eq. ( p5| ) is replaced by the algebraic equation 

= ^ - 1. (103) 

The masses are now measured in units of A^/'^M'^Jm, while N is measured in units of A^/'^M'^Jm? . This implies 
that one must insert a factor of A^/^ in Eq. ( |102| ) to find the total charge on these stars in SI units. 

The curves in Figure 5 are qualitatively similar to those in Figure 3 and the weak field value qmax is the same as in 
the A = case, since the limit imposed by Eq. ( |100D is independent of A. However, the critical charge has decreased 
to a value qc ~ 0.905. This is to be expected, since a factor of a'^A appears in the source term of the scalar wave 
equation ( p5| ) which implies that a large value of A leads to larger inhomogeneities in $. Nevertheless, this effect is 
still fairly minimal. Numerical calculations for other choices of lo show similar behaviour. 

Figures 6 and 7 are derived from the numerical output used to generate Figure 4. Figure 6 is a bifurcation diagram 
for the stars. The first five curves are for solutions with q < qc and each shows a cusp at which the solutions become 
unstable. The remaining three curves, with q > qc, show no cusps and all are stable. The other solution sets we have 
analysed have similar diagrams and all of the solutions with q > qc are stable. However, these objects have such large 
total charges and such small binding energies that their formation may be halted by their own Coulomb repulsion. 

Figure 7 shows the behaviour of the coupling parameter a := —S/M^ for the power law boson stars. Note that 
a here is different from the weak field metric potential defined in Section 0. This coupling parameter measures the 
strength of the coupling between the scalar field and the normal matter (in this case the bosonic matter). In the 
GR limit, a — > 0. As the Figure shows, this quantity is non-negligible even though the theory satisfies the current 
observational constraints. 

The coupling parameter a was first introduced by Damour and Esposito-Farese in the study of neutron star 
equilibrium solutions JlBf . They found that a phenomenon known as spontaneous scalarization (SS) develops in the 
solutions: beyond a certain state of compactness, a ST neutron star develops a non-trivial scalar field configuration 
even when the coupling parameter uj of the ST theory is chosen so that u; — *■ oo far from the star. Their results are 
given in Figure 2 of Rcf. |]l^ and show that, when the baryonic mass Mb of the star is less that some critical value, 
a is negligible, while as Mb is increased beyond this value, a rapidly grows to a maximum before decreasing again as 
Mb is increased to its maximum value. Mb is analogous to Mn for bosons stars and Figure 7 shows that the power 
law boson stars exhibit similar behaviour to neutron stars for large value of M^v but differ considerably when Mn is 
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small. In this case, a increases towards some limiting value as M^v decreases. This limit may be calculated as follows. 



As Mat ^0, (To ^ and from the weak field equations (78) and ([79|) we have 



lim a = - lim f — ) = . (104) 

The difference between these two behaviours may be due to the choice of the coupling parameter aj(<i>). In Ref. [ p^ , 
this coupling has the form 2u; + 3 '--^ l/(log$) which implies that, in the spatially asymptotic region, w ^ oo for the 
choice $oo = 1- For the power law boson stars we have ujoo — 2 /in— 3/2 for = 1, which is finite. 

VII. GRAVITATIONAL EVOLUTION 

When we consider a boson star (or any other compact object) in an evolving cosmological background, we have 
to take into account the influence of the time varying cosmological value of 4> on the structure of the star. To a 
first approximation, and in want of a more rigorous analytical analysis, we assume that $00 increases slowly with 
cosmological time and that the boson star evolves quasi-statically so that, at any point in its evolution, it is described 
by a static equilibrium solution with the appropriate value of ^oo- This approximation is a good one if the free 
fall timescale of the star is much smaller than the characteristic timescale over which G* varies, and is often made 
when considering other compact objects such as white dwarfs. Stars with smaller values of $00 are assumed to exist 
at earlier cosmological times, since in a general ST cosmological model $ tends to increase with time. The star is 
modelled as a sequence of asymptotically flat equilibrium solutions, where each star in the sequence has a value of 
$00 greater than the preceding (earlier time) star. This continues the approach adopted by Torres et al. who 
showed how an increase in affects the physical characteristics of the star. They particularly discuss the evolution 
of a star for which cto remains constant, although they also consider some other cases of more general evolution. 

In this Section we discuss the possible gravitational evolution of both charged and uncharged boson stars of flxed 
particle number N . The assumption that N is constant during the evolution of an asymptotically flat configuration 
is justified since in this case, Eq. ( pO| ) is valid and is a conserved charge. However, in a cosmological setting (in 
which the spacetime is not asymptotically flat) , we cannot in general write a global conservation law for any quantity, 
but instead must approximate a quasi-local conservation equation. This problem is not unique to ST boson stars 
in a cosmological setting: for any metric theory of gravity the formulation of global conservation laws in a general 
spacetime is very difficult. In this analysis we assume that the asymptotic region of the boson star matches smoothly to 
a homogeneous cosmological solution at some large value Xi of the radial coordinate where the boson field amplitude 
a is negligible. We assume that the quasi-local value of the particle number of the star in a cosmological setting 
(evaluated as an integral out to the matching surface x = xi over some space-like hypersurface) is approximately 
equal to N , the particle number the star would have if isolated in an asymptotically flat spacetime. Note that the 
Jordan frame representation of ST gravity that we use in this paper embodies the Einstein Equivalence Principle 
(defined in |jlj]) so that local charge conservation holds and q is constant during the star's evolution. Hence if N is 
constant, this implies that Q remains constant. However, as we shall show below, Mt is not conserved during the 
evolution. Physically, the assumption that N is constant is reasonable since classically, this quantity may be literally 
interpreted as a count of the number of bosons in the star, a procedure that is coordinate independent. The mass, 
on the other hand, is treated as a quasi-local quantity when the star is in a cosmological setting, and there are no 
conservation laws or physical arguments to constrain its behaviour. 

We first examine the evolution of an uncharged boson star of particle number N with A = in BD gravity. We 



make use of the notation of Section IV in which a set of solutions is denoted 5(cro; $00, <?, A), where $00, <1 and A are 
fixed and ctq is a parameter that labels the different elements of S. In the discussion that follows, ctq is an independent 
variable while ai and CT2 are particular values of this variable that label particular solutions within the set. 

Let the initial asymptotic value of the scalar field be <I>i = $00 (ii) at some initial time ti so that the star is a 
member of the set Si{aQ] $1, g = 0, A = 0). We assume that the star is initially stable, which implies that it lies to 
the left of the first maximum in the N{ao) curve |39|. Then the pair (A', $1) uniquely determines the initial central 
amplitude cti of the boson field. The star will have a mass Mi = Mt{<ti, $1, g = 0, A = 0). Here cti is the particular 
value of (To that gives the element of Si with particle number N. At some later time ^2, the asymptotic value of 
the scalar field has evolved to $2 — $00(^2) — k^^i where fc > 1. At this time the star is a member of the set 
iS2((Tq; $2,9 = 0, A = 0) with particle number N. In general, this solution will have a central scalar field amplitude 
(T2 different from ai. From Eqs. (51) and ([55|), the later time solution is generated by a rescaling of the member of 



iSi that has particle number k ^N. Since k > 1, and since for stable stars in Si the particle number curve N{ao) 
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obeys dN/dao > 0, this generating solution has a2 < (Ji- For stable solutions dMr/dag > 0, which implies that the 
mass M2 = MT{a'2,^i,q = 0,A = 0) of the generating solution obeys the relation M2 < Mi. To find the physical 
characteristics of the new solution at time ^2 we must use Eq. (^) to rescale the parameters of the generating solution 
using the rescaling factor k. However, Mt is invariant under this rescaling which implies that the later time boson 
star has a smaller mass than the earlier time solution. Repeating this process for successively higher values of $00 (or, 
equivalently, later cosmological times), one can see that the physical properties of the star during its evolution may be 
calculated from the single initial solution set iSi: as the cosmological time increases, the generating solution is located 
at an ever decreasing value of (Tq and its mass, as well as the mass of the physical solution it generates, decreases. In 
addition, the scalar field becomes increasingly homogeneous and its central value decreases, as viewed from the rest 
frame of the star. Physically, the mass loss may be accounted for by the generation of scalar gravitational radiation, 
which carries energy out from the star, and the weakening of the gravitational coupling strength. As pointed out in 
[Q, within a single solution set the radius of a stable boson star is an increasing function of l/ap (stars approach 
infinite radius in the weak field limit (Tq 0). The radius, like the mass, is invariant under Eq. (|49|), so that this 
quantity also increases with cosmological time. Note that if the star starts out as being stable at time ti, then it will 
remain stable throughout its evolution. 

Figure 8 shows evolution curves for four stable, uncharged boson stars in a; = — 1 BD gravity, where each star has 
a different particle number N. Each star evolves from an initial equilibrium solution with <&oo=l to a configuration 
with $00 = 6. The generating solutions all lie on the q = curves of Figure 2, and all of the stars evolve towards 
lower values of Mt and ap (in other words, the generating solutions approach the weak field limit as the cosmological 
time increases). For the stable part of each set S, the quantity a/1 — which gives a measure of the inverse radius 
of the star, tends towards as ctq ^ (this is true for all of the solutions investigated here). Hence the radius of 
each star increases as it evolves. If we were to continue the evolution beyond $00 = 6, we would find that in the 
limit <I>oo — > 00, the stars tend towards a state of zero mass, zero boson field amplitude and infinite radius (i. e. they 
disperse). Furthermore, since the ratio $0/^00 decreases towards unity as ctq decreases, the scalar field in the interior 
of an evolving star becomes increasingly homogeneous. 

We now compare these results with the evolution curves shown in Figure 1 of Ref. | |4^ . There, it was found that 
the mass of the equilibrium solutions increased with time. This behaviour may be understood if we recall the 
discussion of Section IV. The definition of mass used in |Q differs from the one used here by a factor of $00 • On using 
the definition in the former paper, one finds that the mass increases when we rescale from the generating solution 
(with <I>oo = 1) to the physical solution (with $00 > 1)-|3 The mass we use here is invariant under this rescaling. A 
quantitative example of how these results differ for constant N evolution is given in Table 1. The conclusions in p^ ] 
are seen to be correct once one recognises that their mass is physically different to the definition we adopt here. 

For the sake of completeness, we briefly outline how an uncharged star of constant (Tq evolves. This kind of evolution 
is discussed in Ref. and we reproduce and expand on their results. Since the rescaling factor k increases with 
cosmic time, the generating solutions must be those of ever decreasing values of uo that approach the weak field 
limit as time increases. As is the case for the fixed N evolution, the mass Mt of the constant ao stars decreases 
with time. Here N is no longer conserved and in the physical solution will in general increase with time. This can 
be seen as follows. Under a small change dap in the generating solution we have ap ap — ap + Sap. Since cfq is 
constant in the physical solution, we have fc = 1 — Jcto/cq to first order in Jcto/cq and, since k > 1, this implies that 
6ao < to lowest order. The particle number N changes according to the relation N N — N + Sao{dN/dao), 
where N is the particle number of the generating solution at (Jq and we have written this expression to lowest order 
in daoap. Rescaling to find the corresponding particle numbers of the physical solutions and taking their difference 
we have k'^N — N = Sao{dN/dao — 2N/ao), which is positive wherever the generating solution satisfies the inequality 
{ao/N) dN/dao < 2. One can show numerically that this is true for all stable solutions with q < q^. For any particular 
set of solutions with q > q^, the inequality is satisfied only up to some particular value of (Tq that varies with the 
parameters of the set. Again, as discussed in Section |V, the mass M^ used in Ref. is different from the one 
used here, and it rescales by a factor of fc^ under Eq. p£). In a constant ctq evolution, one can show that will 
be increasing with time if the condition {ao / AI^) dM^ / dap < 2 holds for the generating solution. This is true for 
all solutions with q < qc and for solutions with q > qc whose central boson field amplitude is below some particular 



'^We note here that the choice ^00 = 1 for the generating solution is for convenience only and similar curves could be drawn 
for stars evolving from the distant past (with $00 < 1) through to the present cosmic time (with ~ 1) by making a second 
rescaling of the parameters of the physical solutions. 
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value. 

We next consider the evolution of a charged boson star with A = in BD gravity. The situation here is more 
complex since we must also deal with the charge rescaling in Eq. (^). To start with, we consider a stable boson star 
of particle number N and mass Mi made up of bosons with charge (?, where the asymptotic scalar field has the value 
$1. The solution is a member of the set Si{ao; $i, A = 0). The choice of N and q uniquely determines the central 
boson field amplitude ai of the solution. Again, keeping N fixed we evolve $oo to a new value <I>2 — $00(^2) = 
where fc > 1. The star is now a member of the set 52(cro; <&2, 9, A = 0) with particle number N. However, in this 
case the new solution cannot be generated by rescaling a solution in 5i, since we must also take into account the 
rescaling of q. Instead the generator of the new solution is the member of the set S3{ao;^i,kq,A — 0) with particle 
number k~'^N. Repeating this process for successively higher values of $00 we see that the generating solutions at 
each cosmological time are those of successively higher values q. 

Figure 9 shows the evolution of several charged boson stars in = — 1 BD theory. The stars have charges q = 0.7, 
q = 0.8, q = 0.9 and q = 1.0, and each has = 1 particles. Each of the solutions evolves from an initial value of 
$00 = 1 to a final value of $00 = 2.8. The direction of increasing time in the Figure is from top right to bottom left, 
and each star evolves to a configuration of lower mass and lower central boson field amplitude. Some of the $00 = 1 
solutions used to generate the late time solutions are shown as points in Figure 2, where later time generating solutions 
are those of higher values of q. Their corresponding rescaled (physical) solutions are shown as points in Figure 6. 
Unlike in the uncharged case, charged boson stars cannot evolve indefinitely. At some point in a star's evolution, 
its generating solution has a value of q that exceeds qc and the star will no longer be stable: as the cosmological 
value of G* becomes weaker (ie. <l>oo increases), $ in the interior of the star becomes increasingly more homogeneous 
while its central value also increases. Eventually, the gravitational field inside the star will become too weak for the 
star to remain gravitationally bound and it will be dispersed by its own Coulomb field after some finite time. This 
phenomenon is unique to ST gravity, since it is only in these theories that the gravitational field strength in the star 
may weaken with time. 

The evolution of a star in a more general ST theory is a little more complex since uj will not be invariant under a 
rescaling of $ unless its functional form is changed. For example, given the power law coupling considered above, one 
would also have to rescale the numerical coefficient of 4/3 by a factor of fc~^". Hence, to analyse the evolution of a 
ST star one would have to start with a sequence of generating solutions, each with a different choice of u, much as in 
the case of a charged BD stars in which we needed a sequence of generating solutions, each of a different charge. 



VIII. CONCLUSIONS 



In this work, we have studied the equilibrium configurations of charged boson stars both in Brans-Dicke theory and 
scalar-tensor gravity with a power law coupling, and we have re-examined uncharged boson star solutions. We have 
found that, theoretically, these structures are stable and might exist in the real Universe, provided of course that a 
ST theory correctly describes the gravitational interaction. However, it is likely that only rotating stars will exist in 
an astrophysical setting. While we have made no attempt to construct solutions describing rotating charged stars, the 
existence of rotating stars in the uncharged case leads us to suppose that charged objects of this kind may also exist 
in GR and ST gravity. The use of charged objects to generate astrophysical effects was recently analysed for the case 
of black holes, and this constitutes an additional motivation for the consideration of other charged stellar structures. 
Additionally, the study of the range of possible stellar structures in ST gravity is far from being complete, and a careful 
understanding of the infiuence of an evolving G* in a cosmological background is still missing. Without stating again 
each of the results we have obtained in the previous sections, we may say that the use of simple relativistic objects, 
such as the boson stars explored here, provide a useful setting in which to discuss gravitational issues and to compare 
ST theories with GR. For example, we have found that the properties of both static and evolving charged ST solutions 
in a cosmological background differ considerably from those of GR solutions. By clarifying the mass definitions and 
the properties of the rescaling, we have shown for the first time that gravitational evolution may drive a compact 
object towards a state of decreased energy and central density. The degree to which our conclusions are true for more 
complex and realistic astrophysical objects is unclear. However, it seems likely that these objects will yield observable 
astrophysical signals, like those reported in Refs. (3^,^^ for the case of white dwarfs. Much more work has yet to be 
done, particularly in the field of gravitational memory and gravitational evolution. 
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TABLE I. Here we show some data from an a; = —1, q = BD boson star evolution, in which the particle number is 
conserved. We have chosen N = 0.6. Columns 2, 3 and 4 show the parameters of the generating solutions: central boson field 

amplitude tJo, particle number A'^ and mass (here wc use the Jordan frame ADM mass for ease of comparison). Columns 5 
and 6 show (To and the mass of the physical solution, using the mass definition adopted in the paper by Torres, Schunck 
and Liddlc. Siucc this mass roscalos by a factor , it is larger than in the generating solution and increases with $oo. In 
the present paper, the physical solution has equal ADM mass to the generating solution, since Madm is invariant under the 
rescaling. Thus this mass decreases with cosmic time. 





(To 


N 


Madm 


To 




1.00000 


0.242902 


0.600000 


0.520288 


0.242902 


0.520288 


1.20000 


0.159362 


0.500000 


0.446675 


0.174572 


0.536010 


1.40000 


0.117526 


0.428571 


0.393866 


0.139059 


0.551412 


1.60000 


9.14216E-02 


0.375000 


0.352852 


0.115640 


0.564563 


1.80000 


7.16179E-02 


0.333333 


0.315893 


9.60855E-02 


0.568607 


2.00000 


5.77638E-02 


0.300000 


0.285910 


8.16904E-02 


0.571820 
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0.1 0.2 0.3 0.4 

FIG. 1. Mass curves for the w = 500 Brans Dicke theory with $00 = 1 and A = 0. The curves are parameterised by ao and 
labelled by the boson charge-to-mass ratio q. For each value of q two curves are shown: the Tensor mass Mr (broken curves) 

and the Newtonian mass Mjv (solid curves). Both masses are measured in units of M^j/m. Our choice of $00 implies that Mn 
is numerically equal to the particle number N and is numerically equal to Madm- 
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FIG. 2. Mass curves for the w = — 1 Brans Dicke theory with $00 = 1 and A = 0. Curve labelling and parameterisation are 
the same as in Figure 1. The weak field charge limit is c^ax = 2 and the critical charge is qc = 0.90, above which the mass 
curves diverge. The data points show solutions used to generate the evolution sequences in Figure 9. Four sets of generating 
solutions are shown: q = 1.0 (marked by stars), q = 0.9 (marked by circles), q — 0.8 (marked by triangles) and q = 0.7 (marked 
by squares) . The initial solution for each evolution sequence is at Mn = 1 and the later time generating solutions are those of 
higher values of q and lower values of Mn, Mt and ao. See the text in Section VII for further explanation. 
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FIG. 3. Mass curves for the lo = 1 Brans Dicke theory with $00 = 1 and A = 0. Curve labelling and parameterisation are 
the same as in Figure 1. The weak field charge limit is q^ax = 1-2 and the critical charge is qc = 0.976. 




FIG. 4. Mass curves for the 2a; + 3 = |n$" power law ST theory with n — 4, $00 = 1 and A = 0. Curve labelling and 
parameterisation are the same as in Figure 1. The weak field charge limit is q^ax = 19/16 and the critical charge is qc = 0.986. 
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0.5 1 1.5 2 

FIG. 5. Mass curves for the u! = 1 Brans Dicke theory in the hmit A — > oo with = 1- Curve labelling is the same 
as in Figure 1. The curves are parameterised by the rescaled central density a* and the masses are measured in units of 
Mpi/{mh}^'^). The weak field charge limit is q^ax = 1-2 (the same as in Figure 3) while the critical charge is Qc = 0.905. 
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FIG. 6. Bifurcation diagram for charged boson stars in the n = 4 power law theory. Using catastrophe theory one may show 
that the appearance of cusps signals a change in the stability of the stellar object. Branches that start at the coordinate origin 
are stable. 
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FIG. 7. Scalar coupling a := —S/Mn as a function of rest mass f or ch arged boson stars in the n — 4 Power Law theory. The 
limiting value as Mjv may be obtained analytically using Eq. (|l04). 




FIG. 8. Evolution curves for uncharged boson stars in u — ~1 Brans Dicke theory with A = 0. The curves are labelled by A*' 
and evolve from configurations with "l>oo = 1 to configurations with $oo = 6. The direction of increasing time along the curves 
is towards the origin of the graph. 
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FIG. 9. Evolution curves for charged boson stars in w = — 1 Brans Dicke theory with A = 0. All stars have N = 1 and 
evolve from configurations with $30 = f to configurations with $00 = 2.8. The direction of increasing time along the curves 
is towards the origin of the graph. The symbols (squares, triangles, circles and stars) on each curve are the rescaled $tx3 = 1 
generating solutions shown in Figure 2. 
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